Psychological distress and mental health trajectories during the COVID-19 pandemic in Argentina: a longitudinal study

Psychological-distress increased at the onset of the COVID-19 pandemic in Argentina. Longitudinal studies in developing countries are scarce. Particularly, Argentina had one of the longest lockdowns. Differences in preventive measures against the virus spread between countries may differentially affect the mental health of the populations. Here we aimed to characterize distinct psychological-distress and related-symptoms trajectories associated with the pandemic and explore risk/protective factors. In this longitudinal study, data from 832 Argentineans were collected every 3–5 months, between April 2020–August 2021. Mean psychological-distress levels and related-symptoms tended to increase over time. However, latent-class analysis identified four distinct psychological-distress trajectories. Most individuals had consistently good mental health (Resilient). Two classes showed psychological-distress worsening during the initial phase of the pandemic and recovered at different time points (Fast Recovery; Slow Recovery). The remaining class maintained a mild -level of psychological-distress and began to deteriorate in March 2021 (Deteriorating) continuously. Individuals who are younger, female, have pre-existing psychiatric diagnoses, or have high neuroticism or lower resilience were more likely to experiencing fluctuations in psychological-distress. The mental health trajectory during the pandemic had a complex dynamic. Although most participants remained resilient, a vulnerable group was detected, which deteriorated over time and should be considered by health-services.

In early 2020, the appearance of the novel Coronavirus Disease 2019 (COVID-19) rapidly spread and caused severe health and economic consequences around the globe [1][2][3] . The Covid-19 pandemic and lockdown measures represented a severe and sustained stressor for societies, which raised severe concerns regarding its effects on general well-being and mental health 4,5 . Several cross-sectional studies around the world showed that social isolation and lockdown measures at the initial phase of the pandemic worsened mental health and increased psychological distress, such as higher rates of anxiety, stress-related symptoms, and depressive symptoms [6][7][8][9][10][11] . A recent report estimated a global increase in the prevalence of major depressive disorder of 27.6% and 25.6% for anxiety disorders 12 . Common risk factors for poorer mental health during the early stages of the pandemic were proposed, such as being a woman, being young, having a previous mental health diagnosis, lower education, and a lower socioeconomic status 10,[13][14][15][16][17] . Moreover, trait-measures such as neuroticism, copying-style and resilience, and state-measures such as COVID-19 related fear were found to modulate mental health outcomes 13,[18][19][20] .
As psychological distress may represent an adaptation to environmental threats 21 , it is unknown whether mental health worsening during the initial phase of the COVID-19 pandemic could be transitory or chronic. There is evidence that sustained stressors such as natural disasters may have long-lasting consequences on mental health 22,23 . However, other studies showed that most people were resilient in the long-term, or their mental health improved following an initial deterioration 24,25 . These previous findings are consistent with the idea that stress promotes adaptation and that most individuals are able to cope with negative or threatening experiences. Longitudinal studies performed in developed countries (i.e. Germany, USA, or United Kingdom) within the first six months of the pandemic, reported that the highest levels of anxiety, depression and psychological distress were observed during the early phases of the pandemic and lockdowns [26][27][28][29] . Notably, psychological distress levels tended to decrease or return to pre-pandemic levels as the pandemic unfolded or the restrictive measures were relaxed [26][27][28][30][31][32][33] [but not all studies found this trend [34][35][36][37] ]. Analysis of the trajectory of symptoms over time showed that most individuals were resilient or recovered within the first six months of the pandemic [26][27][28]30 . In contrast, another group of individuals showed sustained poor mental health or deterioration over time, suggesting that they may be at risk in the aftermath 26,30 .
It is thought that the impact of the pandemic may be more significant in developing countries as social inequalities increase, vaccination is delayed and lockdowns are extended 38 . Particularly Argentina had one of the strictest and longest lockdowns in the world. Despite most of the activities being suspended or reduced, the number of new cases grew steadily, and Argentina became one of the countries with the highest rate of infections per inhabitant (Supplementary Material Fig. S1). Previous cross-sectional reports found an initial mental health worsening similar to other countries (i.e. higher rates of psychological distress, anxiety, depression, etc.) and that distress-symptoms were best clustered by severity instead of types 7,13 . No previous study monitored mental health change in time in developing countries or periods longer than 6 months in other locations.
Here, we used a longitudinal sample (n = 832) tracked every 3-5 months since the initial Argentinian lockdown, with the overall aim of analyzing mental health trajectories during 16 months of the pandemic (April 2020-August 2021). More specifically, we aimed to identify distinct psychological distress trajectories over time and characterize associated risk/protective factors.

Method
Participants. We collected data on a convenience sample of 916 Argentine volunteers ranging from 18 to 90 years using an online questionnaire. Participants were recruited using social media, institutional emails, and announcements. Data collection started on April 2020 (Time 1), 11 days after the beginning of mandatory quarantine, continued on July 2020 (Time 2), October 2020 (Time 3), March 2021 (Time 4), and was completed on August 2021 (Time 5). When any of the participants did not complete one or more of the time points, they were excluded from the analysis. In Time 5, the final sample comprised n = 832 participants, representing an attrition rate of 9.2%. Participants did not receive any compensation for their participation. The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008. Furthermore, all procedures were approved by the FLENI ethical committee. Online informed consent was obtained from all subjects.
Instruments. Psychological distress. The 53 item version (range 0-4) of the Brief Symptom Inventory (BSI-53 39 ) was used to assess general mental health along its 9 symptomatic dimensions (Somatization, Anxiety, Phobic Anxiety, Obsession-Compulsion, Interpersonal Sensitivity, Depression, Hostility, Paranoid Ideation, and Psychoticism). Psychological distress was estimated according to the Global Severity Index (GSI) which corresponds to the mean score of each dimension of the BSI-53 and four additional items. The BSI-53 has been used in various psychiatric and natural settings [39][40][41] . It has a 9-factor structure 39,41,42 with robust reliability (α = 0.88). However, this study did not use a representative sample or random population and therefore prevalence levels of GSI or the symptoms dimensions are not presented.
Trait-measures. Big Five Inventory-10 43 was used to assess the big five personality traits: Extroversion, Agreeableness, Openness, Neuroticism, and Conscientiousness. The BFI-10 has a similar structure to the full BFI with acceptable psychometric properties (α = 0.85 44 . BFI-10 employs two items for each trait on a 1-5 Likert scale. Trait-resilience was measured with the 10-item (range 0-4) self-rated Connor-Davidson Resilience Scale 45 . It has a one-factor solution, good reliability (α = 0.91) and validity in non-clinical and clinical samples 46 .
State-measures. COVID-19 related fear (8-items) and Coping skills during the pandemic (5-items), were assessed using two short scales developed in a previous study 13 . Both scales consisted of a 0-4 Likert scale. The COVID-19 related fear scale measures fear of being infected, anxious appraisals about the virus and its potential consequences. On the other hand, the Coping skills during the pandemic scale, assess the ability to cope and tolerate social distancing and lockdowns. In our previous study 13 , we performed a Confirmatory Factor Analysis on the COVID-19 related fear and Coping Skills items and found a two-factor solution with acceptable reliability indices (α = 0.89 and α = 0.79).
Sociodemographic data and covariates. Self-reported sociodemographic data were obtained, including age, gender, occupation, education level, marital status, and income level (Supplementary Material Table S1). Additional covariates were also evaluated: belonging to a known risk group for COVID-19 (yes/no), being an essential service worker (yes/no), having economic concern derived from COVID-19 (range 1-5), the overall number of hygiene measures against COVID-19 (range 1-5), time spent in COVID-19 related information and news (media exposure, range 1-5), importance given to COVID-19 related information and news (media valuation, range 1-5), exercise (yes/no), religiosity (yes/no), tobacco use (yes/no), alcohol consumption (yes/no), being previously exposed to trauma (yes/no) or diagnosed with a neurological or psychiatric disorder (yes/no). Finally, Social Network Size/Strength was measured using the Lubben Social Network Scale 47 , consisting of 12 items on a 7-point scale. This scale has shown robust reliability (α = 0.89).
Statistical analysis. Data analysis was implemented in R, 4.0.5 (R Foundation). When appropriate, categorical and normally distributed variables were analyzed through chi-square tests and ANOVA. Non-normally distributed variables were analyzed with Mann-Whitney-U and Kruskal-Wallis tests.
We first analyzed overall changes in Psychological Distress (GSI) along with Somatization, Anxiety, Phobic Anxiety, Obsession-Compulsion, Interpersonal Sensitivity, Depression, Hostility, Paranoid Ideation, and Psychoticism in Time (April 2020 to August 2021) including different predictors, using hierarchical mixed-effects models. Hierarchical models were implemented using the lme4 package. We used different models that varied in complexity and number of fixed effects for each mental health outcome to evaluate its importance, using participants' ID as a random effect. Model comparison was based on comparing models of different complexity using a Likelihood Ratio Test. When models did not have a significant difference, we then favored less complex models with lower Bayesian Information Criterion (BIC) values. In addition, we reported marginal R 2 and conditional R 2 as effect sizes for the winning models. Marginal R 2 describes the variance explained by the fixed factors while conditional R 2 indicates the proportion of variance explained by the entire model (both the fixed and random factors). Then, we used the Psychological Distress scores (GSI) to construct latent class mixture models (Growth Mixture Models-GMM) and to identify the optimal number of distinct trajectories. GSI was selected to estimate latent classes as it represents an overall index of psychological distress 42 . GMM was implemented in the lcmm package. We used a forward approach, starting from simple models (one-class) to more complex ones (six-class). Moreover, to improve interpretability, age and gender were included as covariates in trajectory analysis. Model fit was assessed using the BIC and entropy levels. After model selection, participants were classified according to their most likely group/trajectory. Finally, to determine potential risk/protective factors associated with each class/trajectory, all covariates and trait/state measures were entered into a mixed-effects multinomial logistic regression. Regression coefficients are presented as odds ratios (ORs) and 95% Cis.

Results
On average, participants were aged 49.8 (SD = 16.4) years old, most of them were women (79.8%), and 24.5% had a previous neuropsychiatric diagnosis (Table S1). Analysis of mean scores over time of Psychological Distress (GSI) and the BSI symptom dimensions revealed to be heterogeneous, as some aspects of mental health improved and others declined (Fig. 1). Model comparison of each symptom dimension revealed that models which included age, gender, and being previously diagnosed with a neuropsychiatric disorder, were those with the better model fit (see Supplementary Material for details). Notably, the inclusion of changes in income, media exposure and valuation, hygiene measures, work changes, and having economic concerns across time, did not improve model fit. In general, individuals who are younger, female, or have pre-existing psychiatric diagnoses  Table S2 for full details).
We then fitted different Growth Mixture Models (GMM) based on the GSI, with one to six latent classes to identify distinct psychological distress trajectories over time. We considered the four-class solution the bestfitting model as it provides lower BIC values than the one-to-three class solution and it is more parsimonious than the five or six-class solution (see Supplementary Material Table S3). The five and six-class solution have the lower BIC values; however, they contain inadmissible sample sizes (< 5% of the entire sample 48 ; Model details of the chosen four-class solution can be found in the Supplementary Material Table S4). Analysis of participant's characteristics showed that age, income level, being previously exposed to a trauma or being previously diagnosed with a neuropsychiatric disorder were significantly different between classes/ trajectories (Table 1).
To determine risk and protective factors associated with class membership, we performed a multinomial logistic regression, using the Resilient class (Class 1) as reference ( Fig. 3; see Supplementary Material Table S5 for details). We found that classes with no consistently good mental health (Fast/Slow Recovery and Deterioration trajectories) were associated with higher scores in trait neuroticism and gender (women). Trait-Resilience and Social-Network size were protective factors against fluctuating (Fast/Slow Recovery classes) or declining (Deteriorating class) psychological distress over time. Moreover, having a previous psychiatric/neurological diagnosis and being previously exposed to trauma were positively associated with Fast Recovery and Slow recovery trajectories. In contrast, adults (30-65 years) and older adults (> 65 years) were associated with lesser odds of being in the Fast Recovery and Deteriorating groups. Sociodemographic variables such as high-income level, were specific predictors of the Deteriorating class. In addition, higher scores of Agreeableness, Conscientiousness, and Openness were associated with a decline of mental health in time (Deteriorating class).
Finally, participants´ most likely latent class membership from the GMM analysis was used to analyze each symptom dimension over time (Fig. 4). That is to say, each of the four class membership was used to analyze the trajectory of the symptom dimension over the 5 time points. As expected, symptom trajectories based on latent classes revealed similar trends as those observed in Psychological Distress trajectories. That is: (1) participants in the Resilient group (Class 1) showed relatively stable low scores over time; (2)   In addition, changes in work, income level, or having economic concern over time did not show a distinct trajectory as a function of class membership (Supplementary Material Fig. S5).

Discussion
In the present study on 832 Argentinian participants, we examined Psychological Distress temporal dynamics across 16 months (April 2020-August 2021) and several symptom dimensions associated with the COVID-19 pandemic. Mental health was revealed to be heterogeneous as fear-related symptoms (COVID-19 fear and Anxiety levels) decreased following the introduction of a strict lockdown, whereas mood symptoms, negative affect, difficulties in social relationship and regulating impulses, tended to increase (Depression, Hostility, Interpersonal Sensitivity, Obsession-Compulsion, and Paranoid-Ideation). However, as the observed effect sizes (R 2 ) are relatively small, these effects should be interpreted cautiously and take into account the wide variability of the pandemic experience between subjects. Furthermore, we identified that being young, women or having a previous neuropsychiatric diagnosis were predictors of more intense psychological distress and related-symptoms over time. Previous longitudinal studies mostly focused on anxiety and depression symptoms and reported an opposite pattern, which indicated that overall mental health improved months after lockdown [26][27][28][29][30][31][32] . These data come from more developed countries (i.e. Germany, England, USA, China) where lockdowns and preventive measures were shorter or more flexible compared to Argentina. Moreover, the economic and social cost between countries is radically different. However, other longitudinal research showed mental health worsening during the first months of the COVID-19 outbreak [34][35][36][37] , suggesting that the relation between the pandemic and Psychological Distress is not uniform but rather a complex phenomenon. Despite finding that mental health deteriorated during the pandemic, the heterogeneous course of symptoms over time may bias the results and obscure the existence of subgroups in the sample. Therefore, trajectory analysis allowed us to analyze the complexity and specificity of the impact of the pandemic among the different individuals. Using latent class analysis, we identified the existence of four distinct symptom trajectories: (1) The Resilient trajectory, which compromised the majority of participants (73.9%), had a consistently good mental health throughout the pandemic as their psychological distress levels and related symptoms were the lowest and stable over time; (2) The Fast Recovery trajectory (10.8%) had a substantial improvement in Psychological distress levels 6 months after the mandatory lockdown (October 2020) where restrictive measures began to relax, and continued to recover until August 2021; (3) The Slow Recovery trajectory (6.7%) displayed a rebound in Psychological distress levels between April 2020-October 2020, when the first wave peaked and showed better mental health following 10 months (August 2021); (4) The Deteriorating trajectory (8.5%) maintained a mild level of Psychological Distress from April 2020 to October 2020 and began to deteriorate in March 2021 until August 2021 where symptoms levels reached their peak. Several circumstances might be related to the improvement  www.nature.com/scientificreports/ of mental health in the Slow/Fast recovery classes over time. For example, relaxation of restrictive measures, the start of the vaccination campaign, outside activity permissions, and school re-openings in some districts.
In addition, mental health improvement over time may indicate a process of adaptation to long stressful events that promote the emergence of new coping skills 25,49 .
Our study described potential risk and protective factors associated with distinct Psychological Distress trajectories. First, individual characteristics such as being young (< 30 years), women, having a smaller socialnetwork size, having a psychiatric/neurological disorder, or being previously exposed to trauma, were associated with an initial increase in Psychological Distress and related symptoms, after the mandatory quarantine (April 2020), followed by a Fast/Slow improvement over time (Fast Recovery and Slow Recovery trajectories). Accordingly, having upper-middle and upper-income levels, were specifically associated with mental health www.nature.com/scientificreports/ worsening over time (Deteriorating trajectory). Second, following previous cross-sectional reports, we provided evidence that trait-characteristic are differentially linked to mental health outcomes 13,18,20 . In this sense, higher levels of neuroticism and lower levels of resilience were associated with higher odds of being in the non-resilient trajectories, and Agreeableness, Conscientiousness, and Openness were specifically related to the Deteriorating trajectory. These results may seem unusual given that high Neuroticism is a strong predictor of worse mental health outcomes, but high Conscientiousness, and Openness, are associated with better outcomes 50 . However, a recent longitudinal study performed in the COVID-19 period found also that high Openness, Agreeableness and Conscientiousness predicted mental health deterioration over time 51 . These results could be understood by considering how the pandemic restrictions affected these behavioral traits. Lockdowns, social distancing, and other economic restrictions, limited the opportunity to engage in new experiences (Openness trait), fulfill altruistic, and social needs (Agreeableness trait), and maintain motivational stability and the fulfillment of responsibilities (Conscientiousness trait). Current findings are in line with other longitudinal reports during the COVID-19 pandemic which found Resilient, Recovered, and Deteriorating trajectories, associated with similar sociodemographic risk factors 26,28 .
The present study has several limitations. Our sample was unrepresentative of the population which may have compromised the observed effect sizes. Most of our respondents were women and participants with middle-to-upper income levels, which may have obscured the relation between economic variables (income and work-related changes over time) and mental health trajectories. In addition, despite the economic efforts made by the Argentinian government, the poverty rate in the first semester of 2021 was 40.9%, which represented approximately a 10% increase compared to pre-pandemic levels (National Institute of Statistics and Census of Argentina-INDEC). Another important limitation is the absence of baseline or pre-pandemic mental health measures that would have facilitated the interpretation of the symptoms trajectories.
To the best of our knowledge, this is the first study to explore the pandemic effects in mental health over a long period of time (> 6 months), and more specifically, in a developing country. With its social and economic inequalities, Argentina had one of the longest and strictest lockdowns in the world, which may have contributed www.nature.com/scientificreports/ to the observed differences between this study and others performed in more developed countries. Overall, we found that mental health trajectory during the COVID-19 pandemic had a complex and heterogeneous dynamic, suggesting that the pandemic experience was different in each individual. At first sight, Psychological Distress levels and related symptoms increased as the pandemic in Argentina unfolded (April 2020-August 2021). However, a careful latent class analysis enabled the characterization of distinct mental health trajectories with their associated risk/protective factors. Finally, this study helped to identify a vulnerable subgroup of individuals with deteriorating mental health, which may need professional assistance. Preventive interventions might be useful to target this group at risk and improve their mental health.

Data availability
The datasets obtained in the current study are available from the corresponding author upon request.